home *** CD-ROM | disk | FTP | other *** search
/ SGI Developer Toolbox 6.1 / SGI Developer Toolbox 6.1 - Disc 4.iso / lib / mathlib / libconv / TRY / simple_diirm1d.f < prev    next >
Encoding:
Text File  |  1994-08-02  |  1006 b   |  38 lines
  1.       subroutine simple_diirm1d(   n_fx, ny, n_gx, f, g, h)
  2.       double precision f(0:(n_fx-1),ny)
  3.       double precision g(0:(n_gx-1))
  4.       double precision h(0:(n_fx+n_gx-2),ny)
  5.       integer n_fx, n_gx, ny
  6. c-----------------------------------------------------------------------------
  7. c
  8. c   This Fortran Subroutine written by
  9. c    Jean-Pierre Panziera
  10. c    Silicon Graphics Inc
  11. c    September 12, 1991
  12. c
  13. c-----------------------------------------------------------------------------
  14.     integer j, ihx, ifx, igx, imax
  15.     double precision zero
  16.     parameter ( zero = 0.)
  17. c-----------------------------------------------------------------------------
  18. c
  19. c   Compute the N 1D convolutions
  20. c
  21. c-----------------------------------------------------------------------------
  22.     do j = 1, ny
  23.         do ihx = 0, n_fx+n_gx-2
  24.         h(ihx,j) = f(ihx,j)
  25.         end do
  26.         do ifx = 0, n_fx-1
  27.         imax = max(n_gx-1,ifx)
  28.         do  igx =  1, imax
  29.             h(ifx,j) = h(ifx,j) - h(ifx-igx,j) * g(igx)
  30.         end do
  31.         end do
  32.     end do
  33.  
  34.       return
  35.       end
  36.  
  37.  
  38.